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Abstract 

Background: Molecular epidemiology of measles virus (MV) is important, not only to measure the success of 
measles vaccination programs but also to monitor the circulation and elimination of the virus worldwide. In this 
study, we compared MV obtained from patients before the 2003 mass vaccination MR campaign and viruses detected 
after 2003 until 2008 in Iran. 

Methods: The nucleoprotein (N) gene of 29 MV strains circulating in Iran between 2002 and 2008 were amplified by 
RT-PCR and subjected to sequence and phylogenetic analysis. 

Results: Molecular characterization of MV studied here revealed that although the outbreaks in Iran were associated 
with MV genotype D4, the isolated viruses clearly belonged to several different lineages. Maximum and minimum 
homology within the 29 Iranian strains in our study was 100% and 94.9% within the carboxyl terminus of the N gene, 
respectively. Using ClustalX program, the alignment of Iranian MV sequences showed nine lineages. 
Conclusion: This study provides the usefulness of MV sequence analysis for the demonstration of local interruption 
of indigenous strain transmission as well as providing a valuable means for monitoring the elimination processes of MV 
control. 
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Introduction 

Measles virus (MV), an enveloped RNA virus 
classified in the family Paramyxoviridae, genus 
Morbillivirus (1), is the most transmissible virus 
known in humans (2). Despite the vaccination pro- 
grams, MV remains a heavy public health bur- 
den worldwide, particularly in developing coun- 
tries (3, 4). Because of the high infectivity of 
measles, 95%-98% vaccine coverage is required 
to prevent virus circulation (2, 5). Serologically 
MV is known as a monotypic virus (6), however 
sequence analysis of the complete H gene and 
the 450 C-terminal nucleotides of the N gene (7, 
8) classifies MV strains into eight clades (A-H) 
dividing into 23 genotypes (9,10). The MV clades 



have been recognized in different parts of the 
world however five of them (Bl, Dl, E, F, and 
Gl) considered inactive since they have not been 
detected in the past 15 yr (11). In the Eastern 
Mediterranean region with established measles 
elimination goals by 2010 the objective has been 
to achieve and maintain interruption of indige- 
nous measles transmission (3). Molecular epide- 
miology of MV is important for either measuring 
the success of measles vaccination programs or 
monitoring circulation and elimination of the virus 
throughout the world (12). Hence, genetic char- 
acterization of wild type MV can help to meas- 
ure transmission pathways and to clarify epide- 
miological links during outbreaks, especially in the 
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countries with an advanced elimination program, 
where most cases are imported (12, 9). 
Measles vaccination started in Iran in 1967. In 
1970 and 2003 the incidence of measles has been 
reported 346/100,000 (13) and 18/100,000 (11644 
measles cases) respectively (14). Despite the re- 
duction in the incidence of disease and the high 
vaccination coverage with routine measles vac- 
cination program, including a two-dose vaccina- 
tion regimen at nine months and 15 yr of age, 
measles cases have been reported in different 
areas in Iran (15). In this context, on December 
2003, Iran conducted a measles/rubella (MR) 
catch-up campaign for individuals 5 to 25 yr olds 
that reached more than 33 milhon people with 
measles vaccine (16). The measles vaccination 
coverage has increased to 98% and the number 
of measles cases has clearly reduced after vac- 
cination program. A post campaign serum-survey 
conducted in 2004 revealed >97.4% of the po- 
pulation aged between 5 and 40 yr had immunity 
to measles (15). Case-based surveillance for mea- 
sles identified three children with laboratory con- 
firmed disease in 2004, 35 in 2005, and 42 in 2006. 
The interruption of indigenous virus circulation 
by mass vaccination campaigns could be dem- 
onstrated by comparing the variability of pre and 
post campaign viruses (5). Importation of a new 
genotype into a country with high immunity to 
measles may lead to outbreaks in the susceptible 
population and genetic analysis will confirm the 
chain of transmission (17, 18). 
In this study, we compared MV obtained from 
patients before the 2003 mass vaccination MR 
campaign and viruses detected after 2003 until 
2008 in Iran. 

Materials and Methods 

This study was approved by the Disease Man- 
agement Center in the Ministry of Health in Iran. 
Patients with rash signed informed consent. Blood 
(serum), throat swabs as well as nasopharyngeal 
specimens, and in some cases, urine that were 
obtained <7 d after onset of a rash from suspected 
cases. All collected samples were transported to 



the National Laboratory of Measles of Iran (NLM) 
in School of Public Health, Tehran University of 
Medical Sciences to process using standard pro- 
cedures. The criteria used for the diagnosis of 
measles and sample collection were as described 
previously by WHO protocols for laboratory di- 
agnosis of measles viral infection (19). MVs were 
detected from throat swabs, nasopharyngeal and 
urine specimens (that were serologically confirmed) 
using B95a and Vero-hSLAM cell lines or directly 
from clinical specimens. RNA was extracted by 
Nucleospin RNA Virus kit according to the 
manufacturer's instructions (M&N, Germany). 
Genotyping of the carboxyl terminus of the nu- 
cleoprotein (N) gene was performed by Reverse 
transcription (RT) and Polymerase Chain Reaction 
(PCR) according to the standard procedures for 
RT-PCR(IO), which is based on a two-step semi- 
nested PCR. Briefly, following reverse transcrip- 
tion, specific primers targeting N gene were de- 
signed to amplify a 599 bp fragment through the 
first round primers (MV60-MV 63.3) using Qiagen 
RT-PCR kit. Second round of PCR were per- 
formed to amplify a 424 fragment using primers 
MV 61 and MV 63 (10). First round PCR prod- 
ucts were sequenced using the ABI prism Big 
Dye Terminator Cycle Sequencing Ready Reac- 
tion Kit (PE Biosystems, Langen, Germany) by 
an Iranian Company (Genefanavaran Co). Geno- 
type identification of the sequences was per- 
formed according to protocols by comparison to 
the WHO reference sequences, representative of 
the 23 established MV genotypes (3, 7, 19). The 
sequencing data was analyzed by Clustal X (20) 
and BioEdit version 7.0.0 DNA analysis software 
(21). Phylogenetic trees were constructed using 
Ttreecon package version 1.3 b (22) using the nu- 
cleotide Kimura-2 parameter and the neighbor- join- 
ing method with bootstrap analysis ( 1 000 replicates). 

Results 

Clinical specimens from patients with acute mea- 
sles infections were collected. Samples belonged 
to 2002-2008, one year before and 5 yr after the 
MR campaign in 2003. A two-step semi-nested 
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RT-PCR was successfully performed and se- 
quences were derived from the most variable re- 
gion of the MV genome encoding the C-terminal 
part of the N-protein of the circulating viruses in 
Iran. Comparison of the analyzed 29 sequences 
with the sequences of the current WHO refer- 
ence strains (7) clearly showed that all Iranian vi- 
ruses in this study were members of genotype D4. 
The derived sequences were registered in Gen- 
Bank (accession numbers: EU234465-EU234487 
and FJ985774- FJ985779). Features of the 29 
sequences are shown in Table 1 . 
Analyzing of the 29 Iranian strains showed the 
maximum and minimum homology within the car- 
boxyl terminus of the N gene was 100% and 94.9%, 
respectively. In the other genotypes (e.g. within clade 
H in China) similar variability was observed (23). 

Table 1: Features of the sequences representative 



It is noteworthy that for phylogenetic analysis we 
used other five Iranian MV strains (lineage 4 and 
5) which were isolated in 2003 before MR cam- 
paign and sequenced in Measles Regional La- 
boratory in Tunisia before measles mass vaccina- 
tion. The result of that study was published by 
Djebbi et al in 2005 (14). 

Using ClustalX program, the alignment of the 34 
adequate sequences showed that the sequences 
from 2002 to 2008 represented in nine lineages. 
According to our findings, five lineages of the 
Iranian strains (lineages 1-5) were found before 
and 4 lineages (lineages 6-9) were detected after 
the MR campaign. The 9 lineages distinction (des- 
ignated 1-9) defined by year of isolation, boot- 
strap values and long branch lengths are shown 
in Table 1 and Fig. 1. 

of 29 sequences obtained during 2002-2008 
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Table 1: continued. 
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Fig. 1: Phylogenetic comparison of 34 Iranian MV strains (C-terminus of N gene) collected from 2002-2008 with other wild- 
type measles genotype D4 detected in other countries. Bootstrap values (1000 replicates) >70% are indicated. The Iranian 
sequences before mass vaccination are shown in bold (green and brown) (Lineages 1-5). The Iranian sequences after mass 
vaccination are shown in bold (blue, purple and red)(Lineages 6-9). NOTE: The Iranian sequences that indicated in bold 
(brown) and italic were sequenced in Measles Regional Center Tunisia before mass vaccination (Lineage 4, 5) (14) 
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Discussion 

The utility of sequence analysis of MV strains 
circulating in a population is a critical means in 
measles, which helps distinguish between endemic 
transmission and importation of the virus (24, 
25). This study investigates the pattern of MV 
genotypes before and after mass vaccination cam- 
paign in Iran. 

The endemic MV strains circulating in Iran were 
not known prior to 1999 and data regarding 
genotype analysis of MV in Iran is limited. 
Briefly, in 1999 during an outbreak in Tehran 
about 30 throat samples from suspected measles 
cases were tested in Centers for Disease Control 
and Prevention (Atlanta-Georgia) , which two of 
them were positive by RT-PCR and MV isolated 
and identified as genotype D4 from only one 
sample (7). Before measles mass campaign in 
2003 in Iran, the NLM in Tehran was (and still 
is) in charge of working on measles suspicious 
cases and virus isolation. From sporadic cases in 
Tehran, five isolated viruses were sent to Tuni- 
sia for genotyping and all of them were again 
identified as genotype D4. 

Following improvement in the ability of surveil- 
lance and epidemiological investigation, we studied 
all IgM positive samples by RT-PCR. Addition- 
ally, nasopharyngeal, throat swabs and urine speci- 
mens were inoculated on B95a and Vero-hSLAM 
cell lines. Attempts also successfully were made to 
isolate MV, as we could isolate MV from 2 na- 
sopharyngeal samples, 1 throat swabs and 1 urine 
specimen. 

The RNA was directly extracted from 25 clinical 
specimens (16 serums and 9 throat specimens) 
and 4 viruses were isolated from cell cultures. 
After MR campaign in Iran, measles surveillance is 
well established and no measles epidemic has 
been reported since 2003. Based on our inves- 
tigations, the phylogenetic analysis clearly showed 
that all of Iranian measles virus strains were in 
genotype D4 and no other obvious patterns were 
observed. We found that within the same province 
(e.g. Kordestan), strains belonging to different line- 
ages could circulate simultaneously in different re- 



gions at different periods. In addition, one given 
lineage could circulate in different provinces at 
the same time. However, before mass vaccination 
no evidence was available for geographic or tem- 
poral restriction, since identical sequences were 
detected in multiple provinces at the same time. 
Nevertheless, there are reports that multiple line- 
ages could co -circulated at the same time (26). 
In our study, high diversity of the detected D4 
strains suggests the role of indigenous viruses be- 
hind several outbreaks in Iran before mass vacci- 
nation campaign. Additionally, it seems that ge- 
notype D4 (with five transmission chains) has been 
responsible for most measles cases in these regions 
before introducing mass vaccination campaign. 
Our results showed that the five Iranian viruses 
(Sistan.IRA/3 1.07/3, Sistan.IRA/3 1.07/2, Sistan. 
IRA/30.07, Siatan. Iran/45.07/1 and Siatan. Iran/ 
47.07/2) as lineage 7 in 2007 and all four Iranian 
viruses (Mazandaran. Iran/10.08, Mazandaran. Iran/ 
17.08, Khorasan. Iran/20.08 and Tehran. Iran/ 
31.08) as lineage 9 in 2008 (after mass vaccina- 
tion campaign) were closely related to other ge- 
notype D4 viruses which were detected in Iran 
as lineage 6 during 2005 and 2006 (Tehran. IRA/ 
25.06 and Golestan. IRA/48.05) (Fig.l). Mean- 
while they were very closely related to wild MV 
which have been reported in the UK in 2006 
(Fig.l) (27, 28). In addition, our results showed 
that Iranian MVs, Sistan. IRA/31.07/4 and Sistan. 
IRA/31.07/1 were more closely related to MV 
from India (Fig.l). 

Genotype D4 is one of the most widely detected 
genotype worldwide with an endemic circulation 
in southern and eastern Africa and in Asian coun- 
tries such as India, Nepal and Pakistan (14). 
In our study, Lineages 1, 4 and 5, the most fre- 
quent before mass vaccination campaign was and 
were present in six provinces throughout our study 
(Fig. 1). 

It has already been shown that the isolates from 
Iran were divided into two lineages (Lineages 4 
and 5) and were quite different from Syrian se- 
quences with 1.3-1.8% nucleotide differences (14). 
The co-circulation of the different lineages within ge- 
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notype D4 in Iran suggests that multiple chains 
of transmission were present in Iran, as such a 
level of divergence within one genotype has al- 
ready been observed in endemic countries, e.g. 
within genotype D4 in the Indian Ocean Islands; 
Mayotte and Seychelles (29), within genotype 
H2 in Vietnam (24), within genotype D8 in Nepal 
(30), within genotype D5 in Cambodia (26) and 
within genotype B3 in Sudan (31). The available 
data from Iran supports the minimum genetic 
heterogeneity within a single chain of transmission, 
regardless of the epidemiologic circumstances. 
Analysis of the predicted amino acid sequences 
of the N proteins showed that the Iranian viruses 
contained several amino acid substitutions that 
have not yet been observed in any of the wild- 
type MV which have previously been sequenced 
(Data not shown). 

As sequence diversity in the C-terminus of the 
N-gene of Iranian MV strains in our study was 
more than 3% before and after the MR cam- 
paign, the D4 genotype detected after the cam- 
paign could have been imported from neighbor- 
ing countries or spread from imported cases 
which were present before the campaign, but re- 
mained undetectable as most of the laboratories 
in provinces lack the necessary equipment and 
supplies for virus isolation. 
It should be mentioned that although the cover- 
age of measles vaccination in Iran is high, it is 
probably not adequate to prevent infection from 
occurring at individual level in some regions and 
creates conditions that favor the rapid transmis- 
sion of a newly introduced genotypes (32). 
It should be noted that the limited registered 
measles virus strains in different databases makes 
the comparison difficult to make sure whether a 
newly detected strain is indigenous in origin or 
imported. 

In conclusion, our study showed that different 
genetic variants of measles virus were circulat- 
ing before and after the vaccination campaign in 
Iran in 2003, although all strains belonged to the 
same genotype (D4). The analysis of genetic 
distances between MV strains before and after 
2003 suggests that indigenous MV transmission 



has been interrupted and the more recent strains 
were imported from other countries. Thus this 
study provides the usefulness of MV sequence 
analysis for the demonstration of local interrup- 
tion of indigenous strain transmission as well as 
providing a valuable means for monitoring the 
elimination processes of virus control. 
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